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1.0 INTRODUCTION 


This User Manual documents the FORTRAN programs RANDOM3 and 
RANDOM4. They are based on fatigue strength reduction, using a probabilistic 
constitutive model. They predict the random lifetime of an engine component to reach a 
given fatigue strength (see Section 2.0, Theoretical Background). 

Included in this Manual are details regarding the theoretical backgrounds of 
RANDOM3 and RANDOM4, input data instructions and sample problems illustrating the 
use of RANDOM3 and RANDOM4 . Appendix A gives information on the physical 
quantities, their symbols, FORTRAN names and both SI and U.S. Customary units. 
Appendix B and C include photocopies of the actual computer printout corresponding to the 
sample problems. Appendices D and E detail the IMSL, Version 10 *, subroutines and 
functions called by RANDOM3 and RANDOM4 and S AS/GRAPH 2 programs that can be 
used to plot both the probability density functions (p.d.f.) and the cumulative distribution 
functions (c.d.f.). 
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2.0 THEORETICAL BACKGROUND 


Fatigue strength data are usually presented as cycles to failure for each of several 
stress amplitudes, the familiar S-N diagram. Results indicate that for lower stress 
amplitudes the cycles (or time) to failure increases. Thus, a power curve fit through the 
data yields a monotonically decreasing curve. In general, this curve is represented as 

S = [N/CT 1/m ' (6) 

where the primitive variables in this equation are as follows: S is the applied constant 
amplitude alternating stress at failure or fatigue strength, N is number of cycles, C' is a 
material parameter that varies from specimen to specimen and m' is a material constant 3 
Equation (6) can be written in terms of "cycles to reach a given fatigue strength" as 


N = CS' m ’ (7) 

Recently another fatigue strength reduction model has been proposed that takes into 
account the effect of temperature as well as other parameters that affect strength. 4 The 
general form of the constitutive relationships for this model is applied to the constituents of 
high temperature composite materials. Specifically, it is applied herein for the case of a 
single material constituent The mechanical property of interest is fatigue strength which is 
expressed in terms of primitive variables, including the general categories of temperature, 
mechanical cycles and mean stress. For these categories, the relationship becomes 
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where S is the applied constant amplitude alternating stress at failure (fatigue strength) at 


current (or operating) temperature, T, mean stress, g, and mechanical cycle, Nm- So is 
fatigue strength at reference temperature, To (usually room temperature), reference mean 

stress (or residual stress), Go, and reference mechanical cycle, Nmo- Also, Tp is the final 
or melting temperature of the material, Sp is the final or tensile strength of the material, and 
Nmf is the final mechanical cycle or lifetime. Empirical parameters, n, m, and q, are 
determined from available experimental data or estimated from anticipated behavior of the 
particular product term. 5 Note that the term containing mechanical cycles is expressed in 
terms of the log of cycles rather than cycles. This formulation is attractive when Nm and 
Nmo are small compared to Nmf. The equation may be solved for Nm, or the "cycles to 
reach a given fatigue strength." The expression is 
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For values typical of a cast nickel base-superalloy subjected to typical loads and 
temperatures, equation (9) indicates increasing life for decreasing temperature, decreasing 
tensile mean stress, and decreasing applied alternating stress. It indicates decreasing life 
for increasing temperature, decreasing compressive mean stress, and increasing applied 
alternating stress. Therefore, equation (9) predicts observed trends in general. 

Probabilistic analysis, via simulation, yields the distribution of the dependent 
random variable, cycles, N. A probability density function (p.d.f.) of cycles is generated 
using the maximum penalized likelihood method for RANDOM3. For RANDOM4, a 
p.d.f. of cycles is generated using the maximum entropy method. Maximum entropy uses 
Jaynes’ principle which says that "the minimally prejudiced distribution is that which 
maximizes the entropy subjected to the constraints supplied by the given information." 6 


3.0 INPUT DATA 


Data input for RANDOM3 and RANDOM4 is user friendly and easy to manipulate 
(see, for example, the file entided NORMAL.INP, in Section 4.0). The first twelve lines 
of input have die same format, 2E12.4 and the last two lines differ. The last two lines of 
input have the formats I3,2X,I3,2X,2E12.4,2X,I3 and 13, respectively. A brief, line by 
line description is given along with an example for each line (NOTE: the ruler is to aid the 
user in formatting and is not a part of the input). A table listing the physical quantities, 
their units and symbols is given in Appendix A. 

1. Random Number Generator Seed, ISEED, and Sample Size, NTOT 

EXAMPLE: 

123456789012345678901234567890 
1 40 


2. Ultimate Tensile Strength, SF 
EXAMPLE: 

123456789012345678901234567890 

900.0000 45.0000 


3. Log of Final Cycle, NMF 
EXAMPLE: 

123456789012345678901234567890 

8.0000 0.8000 

4. Reference Fatigue Strength, SO 
EXAMPLE: 

123456789012345678901234567890 

500.0000 25.0000 


5. Log of Reference Cycle, NMO 
EXAMPLE: 

123456789012345678901234567890 

7.0000 0.7000 
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6. Current Fatigue Strength, S 
EXAMPLE: 

123456789012345678901234567890 

250.0000 12.0000 


7. Residual Compressive Stress, SIGO 
EXAMPLE: 

123456789012345678901234567890 
20.0000 1.0000 


8. Current Mean Stress, SIG 
EXAMPLE: 

123456789012345678901234567890 

150.0000 7.5000 


9. Temperature Exponent, XXN, Stress Exponent, XXM, and Cycle Exponent, XXQ 
EXAMPLE: 

123456789012345678901234567890 
0.5000 0.0150 


10. Melting Temperature, TF 
EXAMPLE: 

123456789012345678901234567890 

1500.0000 75.0000 


1 1. Reference Temperature, TO 
EXAMPLE: 

123456789012345678901234567890 

20.0000 0.6000 
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12. Current Temperature, T 
EXAMPLE: 

123456789012345678901234567890 
850.0000 25.0000 


13. The DESPL"! parameters are NODE, INIT, ALPHA, EPS, and MAXIT and are 
entered in that order as follows: 

EXAMPLE: 

1234567890123456789012345678901234567890 
21 0 20.0000 1.0E-05 30 

14. The DESPL parameter, IOPT, is entered as follows: 

EXAMPLE: 

1234567890 

2 
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4.0 SAMPLE PROBLEMS FOR RANDOM3 AND RANDOM4 


The objective of these programs is to predict the random lifetime to reach a given 
fatigue strength for an engine component. The theory is based on fatigue strength 
reduction, using a probabilistic constitutive model. The only difference between 
RANDOM3 and RANDOM4 is the method used to generate p.d.f. estimates. RANDOM3 
uses maximum penalized likelihood, while RANDOM4 uses maximum entropy (see 
Section 2.0, Theoretical Background). RANDOM3 and RANDOM4 input parameters are 
given in Table A2.1. 


TABLE A2.1 RANDOM3 and RANDOM4 input (SI units) 


FORTRAN Distribution Mean Standard Deviation 


Name Type (Value) (% of Mean) 


SF 

normal 

900.0 

45.0 

(3%) 

NMF 

lognormal 

8.0 

0.8 

(10%) 

SO 

lognormal 

500.0 

25.0 

(5%) 

NMO 

lognormal 

7.0 

0.7 

(10%) 

S 

lognormal 

250.0 

12.5 

(5%) 

SIGO 

lognormal 

-20.0 

-1.0 

(1%) 

SIG 

lognormal 

150.0 

7.5 

(5%) 

XXN 

normal 

0.5 

0.015 

(0.3%) 

XXM 

normal 

0.5 

0.015 

(0.3%) 

XXQ 

normal 

0.5 

0.015 

(0.3%) 

TF 

normal 

1500.0 

45.0 

(3%) 

TO 

normal 

20.0 

0.6 

(3%) 

T 

normal 

850.0 

25.5 

(3%) 
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The input is entered in the following format in a file 
entitled NORMAL.INP. 


1234567890123456789012345678901234567890 


1 

US 

900.0000 

45.0000 

8.0000 

0.8000 

500.0000 

25.0000 

7.0000 

0.7000 

250.0000 

12.5000 

20.0000 

1.0000 

150.0000 

7.5000 

0.5000 

0.0150 

1500.0000 

75.0000 

20.0000 

0.6000 

850.0000 

25.5000 


21 0 20.00 1.0E-05 30 

2 


Execution of RANDOM3 and RANDOM4 (source code entitled NR3.FOR and 
NR4.FOR, respectively) produces files entitled RANDM33 and RANDM44. These give 
intermediate results (see Appendices B and C). Execution also produces plotfiles entitled 
PLOT1 and PLOT2 (see Appendices B and C). These files are used to plot the X and Y 
axes of the probability density function (p.d.f.) and the cumulative distribution function 
(c.d.f.), respectively, generated by RANDOMS and RANDOM4. The plots are drawn 
from die plotfiles by the SAS/GRAPH graphing program (see Appendix D). These plots 
for the sample problem are shown Figures 1, 2, 3, and 4. This same sample problem has 
been reported in Boyce and Chamis. 7 There, however, it utilized U.S. Customary units 
and older versions of RANDOM3 and RANDOM4 (using IMSL Version 9.2 subroutines) 
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Fig. A2.1 p.d.f. of log of mechanical cycles for fatigue strength reduction model, 
using maximum penalized likelihood method of p.d.f. generation. 



Fig. A2.2 c.d.f. of log of mechanical cycles for fatigue strength reduction model, 
using maximum penalized likelihood method of p.d.f. generation. 
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Fig. A2.3 p.d.f. of log of mechanical cycles for fatigue strength reduction model, 
using maximum entropy method of p.d.f. generation. 



Fig. A2.4 c.d.f. of log of mechanical cycles for fatigue strength reduction model, 
using maximum entropy method of p.d.f. generation. 
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6.0 APPENDIX A 


PHYSICAL QUANTITIES, SYMBOLS, AND UNITS 


The physical quantities, their symbols and units for the fatigue crack growth 
model are given in the following table. 


Table A2.2 Physical quantities, symbols, and units for fatigue 
crack growth model for RANDOM3 and RANDOM4. 


Physical 

Theory 

FORTRAN 

Units 


ESI 

Name 

SI U.S 

Ultimate Tensile 
Strength 

SF 

SF 

MPa ksi 

Final Cycle 
(lifetime) 

n mf 

NMF 

dimensionless 

Reference Fatigue 
Strength 

SO 

SO 

MPa ksi 

Reference Cycles 

n mo 

NMO 

dimensionless 

Current Fatigue 
Strengths 

S 

S 

MPa ksi 

Residual Compressive 
Stress 

o 0 

SIGO 

MPa ksi 

Current Mean Stress 

c 

SIG 

MPa ksi 


n 

XXN 

dimensionless 

Empirical Material 
Parameters 

m 

XXM 

dimensionless 


q 

XXQ 

dimensionless 

Melting Temperature 

TF 

TF 

u* 

o 

u 

o 

Reference 

Temperature 

TO 

TO 

o 

U 

o 

Current Temperature 

T 

T 

°c °F 
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7.0 APPENDIX B 


RANDOM3 SAMPLE PROBLEM: SOURCE, INPUT AND OUTPUT FILES 
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9.0 APPENDIX D 


IMSL SUBROUTINE CALLS FROM RANDOM3 AND RANDOM4 


RANDOM3 

1. RNSET - Initializes a random seed for use in the IMSL 

random number generators. 

2. RNNOR - Generates pseudorandom numbers from a standard 

normal distribution using an inverse CDF 
method. 

3. RNLNL - Generates pseudorandom numbers from a 

lognormal distribution. 

4. DESPL - Performs nonparametric probability density 

function estimation by die penalized 
likelihood method. 

5. GCDF - Evaluates a general continuous cumulative 

distribution function given the ordinates 
of the density. 


RANDOM4 

1. RNSET - Initializes a random seed for use in the IMSL 

random number generators. 

2. RNNOR - Generates pseudorandom numbers from a standard 

normal distribution using an inverse CDF 
method. 

3 . RNLNL - Generates pseudorandom numbers from a 

lognormal distribution. 
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10.0 APPENDIX E 


SAMPLE SAS/GRAPH PROGRAM FOR RANDOM3 AND RANDOM4 


data a; 

INFILE TLOT1.CPR' FIRSTOBS=2;input x y; 

GOPTIONS DEVICE=HP7470; 
proc gplot; 

axisl label=(h=l f=simplex LOG OF CYCLES') 
value=(h=l f=simplex); 
axis2 value=(h=l f=simplex) label=none; 
plot y*x / haxis=axisl vaxis=axis2; 

TITLE H=1 A =90 F=SIMPLEX 'PROBABILITY DENSITY FUNCTION'; 
symbol i=spline v=square; 
dataB; 

INFILE 'PLOT2.CPR' FIRSTOBS=2;input x y; 
proc gplot; 

axisl label=(h=l f=simplex ’LOG OF CYCLES') 
value=(h=l f=simplex); 
axis2 value=(h=l f=simplex) label=none; 
plot y*x / haxis=axisl vaxis=axis2; 

TITLE H=1 A=90 F=SIMPLEX 'CUMULATIVE DISTRIBUTION FUNCTION'; 
symbol i=spline v=square; 
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